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Exclusive proton-antiproton annihilation into two photons at large s (~ 10 GeV 2 ) and \t\, \u\ ~ s 
can be described by a generalized parton picture analogous to the "soft mechanism" in wide-angle 
real Compton scattering. The two photons are emitted in the annihilation of a single fast quark 
and antiquark. The matrix element describing the transition of the pp system to a qq pair can be 
related to the timelike proton elastic form factors as well as to the quark/antiquark distributions 
measured in inclusive deep-inelastic scattering. The reaction could be studied with the proposed 
1.5 — 15 GeV high-luminosity antiproton storage ring (HESR) at GSI. 



Compton scattering — both real and virtual — is one 
of the main sources of information on the structure of the 
nucleon. In particular, inclusive deep -inelastic scattering 
can be viewed as a measurement of the imaginary part 
of the forward virtual Compton amplitude, parametrized 
by the quark and antiquark distributions in the nucleon. 
More recently, exclusive processes have been considered, 
which can be described in terms of generalized parton 
distributions, namely deeply-virtual (DVCS) and wide- 
angle real (WACS) Compton scattering. It has been ar- 
gued that at values of s ~ few GeV 2 and \t\,\u\ ~ s 
WACS is dominated by the "soft mechanism" ^] . The 
Compton scattering occurs off a single quark or antiquark 
in the nucleon. Its emission and absorption by the nu- 
cleon is described by double distributions which can be 
related to the usual quark/antiquark distributions as well 
as to the elastic form factors of the nucleon jD]. This 
approach describes well the existing data ||, including 
the spin asymmetry of the cross section measured re- 
cently in the JLAB Hall A experiment The hard 
scattering mechanism, in which the struck quark rescat- 
ters via gluon exchanges of virtuality ~> t, is relevant only 
at asymptotically large s and t and cannot account for 
the measured cross section at JLAB energies ]5J . 

The proposed high-luminosity 1.5 — 15 GeV antiproton 
storage ring (HESR) at GSI [|| would offer an opportu- 
nity to study the Compton process in the crossed channel, 
namely exclusive proton-antiproton annihilation into two 
photons, pp — > 77. In this letter we argue that this pro- 
cess at large s and |u| ~ s can also be described in 
a generalized parton picture. The two photons are pre- 
dominantly emitted in the annihilation of a single "fast" 
quark and antiquark originating from the proton and an- 
tiproton. The new double distributions, describing the 
transition of the pp system to a qq pair, can be related to 
the timelike nucleon form factors; by crossing symmetry 
they are also connected with the usual quark/antiquark 
distributions in the nucleon. A similar approach based on 
light-cone wave functions has been used for the processes 
77 — > 7T7T jf| and 77 — ► BB . Exclusive pp annihilation 
was also studied in the diquark model H . 



The annihilation process p{p\) +p(p2) ~* 7(91) +7(92) 
is characterized by the invariants s = 2{p\p2) > (we 
neglect the proton mass), t = — s sin 2 (#/2) < 0, and 
u = — s cos 2 (9/2) < 0, where 9 is the scattering an- 
gle in the center-of-mass frame. We consider the region 
where s is much larger than typical light hadron masses, 
s ~ 10 GeV 2 , and \t\, \u\ ~ s, which requires that 6 be 
sufficiently far from and 7r (wide-angle scattering). 
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FIG. 1: The "handbag" contribution to pp —> 77 annihilation. 

In QED, the annihilation process e + e~ — > 77 in- 
volves the t (or u) channel exchange of a virtual elec- 
tron/positron with spacelike four-momentum. In pp — ► 
77 annihilation in QCD, the exchanged system contains 
at least three quarks. At large momentum transfer such 
an exchange is strongly suppressed by the nucleon wave 
functions. In this situation, the most efficient way of ac- 
complishing a large momentum transfer is the "handbag" 
diagram, Fig. [ij The amplitude consists of two parts. In 
the first part, the proton-antiproton system converts to 
a quark-antiquark pair by exchanging a virtual qq ( "di- 
quark" ) system. The nucleon wave functions restrict the 
spacelike virtuality of this system to values of the order 
of a hadronic scale, A 2 , related to the size of the nucleon. 
In the second part, the quark-antiquark pair annihilates 
via exchange of a highly virtual quark/antiquark, much 
like in QED e + e~ — > 77 annihilation. Neglecting trans- 
verse momenta, we expand the momenta of the quark 
and antiquark as 

kl = (1 -X!)p! +X2P2, (1) 

h = X1P1 + (1 - x 2 )p2- (2) 
The variables x\ and X2 obey the "partonic" limits < 



2 



xi.2 < 1. The spacelike virtuality of the exchanged "di- 
quark" system is —x\Xis. The restriction |a^x^2s| < A 2 
implies that for s>A 2 the value of either x\ or x 2 must 
be small, of the order A 2 /s. This means that either the 
annihilating quark or the antiquark are "fast" , carrying 
the major part (~ 1) of the proton or antiproton mo- 
mentum. This in turn implies that the spacelike virtuali- 
ties of the quark propagators connecting the photon ver- 
tices in the graphs of Fig. |l|, [— x\Xis + (1 — x\ — X2)t] for 
the direct and [— X1X2S + (1 — x\ — X2)u] for the crossed 
graph, are large in the kinematics of wide-angle scatter- 
ing, |t|, lit | ~ s. Thus, the parton picture of pp — ► 77 ex- 
clusive annihilation emerging from the "handbag" graph 
is self-consistent. The inclusion of transverse momenta 
would change this picture only quantitatively. 

The so-called hard scattering mechanism would require 
the rescattering of the active partons through gluon ex- 
changes of virtuality t. In the case of WACS such contri- 
butions, which arc theoretically dominant in the asymp- 
totic region |t| — * 00, are negligible compared to the sub- 
asymptotic "handbag" contribution at all experimentally 
relevant values (f ~ few GeV 2 ), the main reason being 
that the hard contributions are numerically suppressed 
by a factor of (a s /ir) 2 w 1/100, see Ref. [Q. The same 
can be expected for the annihilation channel. 

The crucial ingredient in the calculation of the 
hadronic amplitude is the matrix element describing the 
amplitude for the conversion of the pp system to a qq 
pair. Following the treatment of WACS, we parametrize 
this matrix element by "double distributions" |l], |^, (To] , 
in which the quark and antiquark momenta are mea- 
sured in terms of the proton and antiproton momenta 
in a symmetric fashion, using the average and difference, 
P = (pi +P2)/2 and r =pi - p 2 : 



ki = (1 + a)p + xr/2, 
fc 2 = (1 — a)p — xr/2. 



(3) 
(4) 



The new variables x and a are related to the old ones 
by x — 1 — x\ — X2, a — X2 — x\. We parametrize the 
annihilation-type matrix elements of the quark bilinear 
vector operators as (a denotes the quark flavor) 

(0| Vv(-z/2) 7 ^ a (z/2) |pi,Ai;p 2 A 2 ) 

= v la u J J dxda F a (x,a;s) e- iaipz) - lx{rz)/2 
M+M<i 

+ .... (5) 

The matrix element of the axial vector operator (7^75) 
is parametrized similarly with a double distribution 
G a (x,a;s). Here z is a light-like distance, and u = 
u(pi,Xi) and v = v(p 2 ,A 2 ) are the proton and antipro- 
ton spinors. The double distribution F a depends on the 
spectral parameters x and a, as well as on s. Not shown 
in Eq.(^) for brevity are components of the matrix ele- 
ment of the type e CTQ( 3 7 W7 7 75W, which play a role in main- 



taining electromagnetic gauge invariance of the "hand- 
bag" amplitude, similar to the kinematical twist-3 terms 
in deeply-virtual Compton scattering J[l] . For simplic- 
ity we neglect components of the matrix element corre- 
sponding to the Pauli and pseudoscalar form factors; they 
should be included in a more complete treatment. 

The new double distributions are related to the time- 
like elastic form factors of the proton by 



^~ /£a J J ^ X ^ ^ a ^ Xl a ' S ^ 
a M+M<i 



(6) 



where F\(s) is the Dirac vector form factor. A similar 
relation holds for the function G a {x,a; s) and the axial 
form factor. Furthermore, in the limit s — > one can use 
crossing symmetry to relate the annihilation-type matrix 
element (||) to the corresponding scattering-type matrix 
element at t = 0, parametrized by the usual double dis- 
tributions of WACS. In particular, this implies 



da F a (x, a; s = 0) 



fa(x), 



(J) 
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where f a (x) is the usual unpolarized parton density in 
the proton; in conventional notation f a (x) = 9(x)q a (x) — 
6{—x)q a {—x). In a similar way the double distribution 
G a {x 1 a\ s) reduces to the polarized parton density. 

To construct an explicit model for the double distribu- 
tions we make the ansatz 

F a {x,a\s) = j a (x) h a (x,a) S a (x,a;s). (8) 

In the parton density, f a {x), we take into account only 
the valence quark contribution, i.e., we assume it to be of 
the form f a (x) = 8(x) q a (x); sea quarks could readily be 
included. The factor h a (x, a) is a profile function, whose 
integral over a from — 1 + \x\ to 1 — \x\ is normalized to 
unity, cf. Eq. (]?]). Finally, S a (x, a; s) is a cutoff function 
accounting for the s-dependence of the double distribu- 
tion, defined such that S a (x, a; s = 0) = 1. We model it 
as (assuming x > 0) 



SaixT^s) 



exp 



{{l-x) 2 -a 2 ]s 



or, in terms of the original variables x\ and x 2 



exp 



X1X2S 



(Xl + X2){1 - Xl - X 2 )A 2 



(9) 



(10) 



Here A 2 is a parameter of dimension mass squared. This 
cutoff suppresses contributions in which the absolute 
value of the virtuality of the exchanged "diquark" sys- 
tem, |xix 2 s|, is large, and also configurations in which 
xi + X2 tends to 1 (x tends to zero). For the double 
distributions parametrizing scattering-type matrix ele- 
ments in WACS, this cutoff function can be thought of 
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as the result of an overlap integral of light-cone wave 
functions jll ||] . Such an interpretation is not possible in 
the annihilation channel. Here, the cutoff function should 
simply be regarded as an effective manner to restrict the 
virtuality of the exchanged "spectator" system. 

For an estimate of the pp — > 77 amplitude we use the 
GRV94 LO parametrization of the valence quark den- 
sities at a normalization scale of 1 GeV 2 [G3 . We ne- 
glect small contributions from strange quarks. The pro- 
file function we take to be of the form h a (x,a) = S(a); 
extended profiles have been suggested in Refs.j9|, [ToJ . 
For simplicity we choose the cutoff function S a (x, a;s) 
to be the same for both u and d flavors. The parame- 
ter A 2 = A 2 = A 2 we determine by fitting the integral 
over the model double distribution [see Eq.(^)] to the 
data for the proton vector form factor. Fitting to the 
form factor data in the spacelike domain |]l3f we obtain 
A 2 = 0.7 GeV 2 . For large s the real part of the time- 
like form factor is approximately two times as large as 
that of the spacelike form factor for corresponding large 
t, see Ref.Q for a discussion. In order to account for 
this effect we multiply the timelike form factors entering 
the double distribution ansatz by a factor of 2. We stress 
that this should be regarded as a purely phenomenolog- 
ical improvement. The axial vector double distribution 
G a (x,a;s) is modeled analogously, using the GRSV95 
LO parametrization for the polarized valence quark den- 
sities For simplicity we use here the same cutoff 
function as in the vector double distribution F a (x, a; s). 

With our model for the double distributions we can 
compute the pp — > 77 amplitude from the "handbag" 
graphs of Fig. ff| As the process is dominated by con- 
figurations with small x\ and x^ ("fast" quark and an- 
tiquark), we can simplify the qq — > 77 amplitude by ne- 
glecting contributions of order X1X2] keeping them would 
result in A 2 /s-suppressed contributions to the pp — > 77 
amplitude which are beyond the accuracy of our approxi- 
mation. In this way the annihilating quark and antiquark 
are effectively put on mass shell, which makes it straight- 
forward to maintain transversality of the amplitude (e.m. 
gauge invariance). The result for the helicity-averaged 
differential cross section is 



da 
dcos ( 



27rag m R 2 v {s) cos 2 9 + R\{s) 
s sin 2 9 



(11) 



The information about the structure of the proton is con- 
tained in generalized form factors 



R v (s) 



dx da 

x| + |ck|<1 



F a (x,a; s) 



(12) 



and Ra(s) defined by the corresponding integral over 
the double distribution G a (x,a;s). For Ry = Ra = 1 
Eq.(O) would reduce to the Klein-Nishina formula for 
the s~ — > 77 cross section in QED. Note that our par- 
ton picture is applicable only to wide-angle scattering, so 
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FIG. 2: The squared form factors Ry(s) (solid line) and 
R A ( S ) (dashed line ), a s calculated from the double distribu- 
tion model, cf. Eq.(|12j). 



the divergence of the expression in Eq. (|1 1[) in the limit 
9 — ► or 7r should be regarded as unphysical. 

Eq. ([n]) shows that the contribution of the vector oper- 
ator to the annihilation cross section, Ry, is suppressed 
at large scattering angle relative to that of the axial vec- 
tor operator, Ra, by a factor of cos 2 9. This is differ- 
ent from WACS, where the contribution from the axial- 
vector form factor is suppressed 0. 

Compared to the elastic form factor, Eq.(|J), the inte- 
grand in Eq.(|l2l) contains an additional factor 1/x, which 
is the "remnant" of the quark propagator in the qq — > 77 
amplitude. Note that the integral nevertheless converges 
at small the cutoff function S a (x,a;s) forces the 

double distribution to vanish for x — > 0, cf. Eqs.(||) and 
(^). More generally, the properties of the double distri- 
butions ensure that the integrals are dominated by large 
values of x. The numerical results for the squared form 
factors Rv{s) and Ra{s) are shown in Fig. |[ 

The results for the form factors Rv(s) and Ra(s) turn 
out to be insensitive to the precise value of the normaliza- 
tion scale of the parton densities in the double distribu- 
tion model, Eq.(||). This happens due to the correlation 
of the A parameter in the cutoff function with the nor- 
malization scale of the parton densities |l6[| . Changing 
the normalization scale one must change A such as to re- 
fit the form factor data, and the two changes compensate 
each other in the integral Eq. ( |l2"| ) . This fact is important 
for the consistency of our approach. 

In our simple model the timelike double distributions 
are described by real functions. As a result, the pp — » 77 
amplitude is also real (the virtuality of the quark propa- 
gators in the "handbag" graphs of Fig. [j] is always space- 
like). A more refined treatment should include also the 
"intrinsic" imaginary part of the double distributions, 
which is related to the imaginary part of the timelike 
proton form factor by Eq.(^|). 
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FIG. 3: The pp — + 77 cross section integrated over the range 
45° < 6 < 135°, as a function of s. 



(Pauli and pseudoscalar form factor), and transverse mo- 
menta of the partons. Also, the picture proposed here 
can be applied to polarization observables, as well as to 
baryon-antibaryon production in 77 reactions. 
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It is interesting to estimate the counting rate for pp — > 
77 annihilation expected for the proposed 1.5 — 15GeV 
antiproton storage ring (HESR) at GSI [gj. With a 
solid target the luminosity could be as high as L = 
2 x 10 cm s . Since our parton picture applies 
only in the kinematical region where \u\ ~ s, we in- 
tegrate the differential cross section ( |TT| ) over a range 
6 < 8 < 7r — 6, with S > 0, excluding small scattering an- 
gles. For s = lOGeV 2 , a range of 45° < 9 < 135° corre- 
sponds to \t\, \u\ > 1.5 GeV 2 , which seems to be a reason- 
able boundary in view of the experience with wide-angle 
Compton scattering |Q . The cross section integrated over 
this fixed angular range (divided by 2 to account for the 
identical particles in the final state) is shown in Fig. |^ as 
a function of s. At s = 10 GeV 2 our model predicts an 
integrated cross section of 0.25 x 10 _9 fm 2 , correspond- 
ing to a counting rate of 0.5 x 10 _3 sec _1 , that is O(10 3 ) 
events per month. The process should thus be measur- 
able with reasonable statistics. 

To summarize, we have outlined a generalized parton 
picture of exclusive pp — > 77 annihilation, based on the 
"handbag" graph and certain assumptions about the soft 
matrix element describing the conversion of the pp sys- 
tem to a qq pair. The approximations made in the present 
treatment can be refined in many ways, most notably by 
including other components of the soft matrix element 
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